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We show that a vortex in a chiral p-wave superconductor, which has the px + ipy-wave pairing 
state and breaks (7(1), parity and time reversal symmetry simultaneously, has fractional charge — ^ 

and fractional angular momentum — {n\ vorticity). This suggests that the vortex could be anyon 
and could obey fractional statistics. Electromagnetic property of the vortex is also discussed and 
we find that an electric field is induced near the vortex core. 



I. INTRODUCTION 

Recently, the superconjductivity of Sr2Ru04 wa.s dis- 
covered by Maeno et. aL.u This superconductivity seems 
to exhibit chiral p-wave order, in which the symme|t|ry of 
the wave function of Cooper pair is px + ipy-wavea and 
parity (P) and time reversal symmetry (T) are violated. 
In this paper, we show that a vortex in chiral p-wave su- 
perconductor has fractional charge and fractional angular 
momentum as the consequence of P- and T-violation. 

Sr2RuOi has a layered perovskite crystal structure 
without cupper and has low transition temperature 
Tc ~ l.b[K]. Sr2RuCL js a strongly correlated 2- 
dimensional Fermi liquidB^Q and is partially similar to 
•^He which shows p-wave superfluidity.El The evidence 
for unconventional (i.e., non s-wave) superconductivity 
was shown by the absence of|-a Habel-Slichter peak in 
1/TiT of NQR-measurejjpentaD and Tc-suppression by 
non-magnetic impurities.El Rice and Sigrist proposed that 
pairing symmetry of orbital part is chiral p-wave,B which, 
has the same orbital symmetry of superfluid '^He-A.El 
The proposal is consistent with the discovery of inter- 
nal magneticJield in the superconducting phase by /iSR 
measpjemenlO and the experiment of '^^O-NMR Knight 
shift .t3 In the chiral p-wave superconductor, Cooper pair 
has orbital angular momentum Iz — I, where z is per- 
pendicular to the superconducting plane, i.e., the same 
direction as the c axis of the crystal. Therefore, local 
U{1) gauge symmetry, P and T are spontaneously bro- 
ken by the same order parameter. 

The above situation is analogous to the quantum Hall 
system (QHS), which is the 2-dimcnsional electron sys- 
tem in an external magnetic ficld.Ell In QHS the external 
magnetic field violates P and T, but the explicit local 
U{1) gauge symmetry is preserved. It was shown that 
the Chern-Simons term is induced in the effective actiojo. 
of the gauge field when Fermion is integrated out.Ej'tj 
The Chern-Simons term is a P- and T-odd bi-linear form 
of the gauge fields and has one derivative. Therefore, 
this term plays important roles in low energy and long 
distance physics and causes P- and T-violating electro- 
magnetic phenomena. The coefficient of the induced 
Chern-Simons term is quantized exactly as (integer) x 
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U{1) gauge invariance guarantees the exact quan- 
tization. Since the coefhcient of the term becomes the 



Hall conductance, the integer quaritum Hall effect was 
explained from gauge invariance .EJc3 

It was shown that the Chern-Simons term is also in- 
duced in the Gipabyrg-Landau action of chiral p-wave 
superconductors.tatZI In this case, the Chern-Simons 
term does not have exactly quantized coefficient because 
of the spontaneous breaking of local U{1) gauge symme- 
try. The coefficieiit, is approximately equal to the fine 
structure constant Jlj 

In this paper, we discuss a vortex in chiral p- 
wave superconductor starting from the Ginzburg-Landau 
Lagrangianllj with the Chern-Simons term. Such a vor- 
tex could obey unusual statistics by the existence of 
the Chern-Simons term. The above possibility has been 
pointed jXiut in Chern-Simons- Higgs model by several 
authors.E3 In this model, vortex solutions have been 
found and they have fractional angular momentum and 
fractional charge. The Chern-Simons vortex is a can- 
didate of anyon that obeys the fractional statistics due 
to fractional angular momentum.E3 In the present pa- 
per we find a vortex solution in chiral p-wave supercon- 
ductor and show that it has fractional charge — ^ and 
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fractional angular momentum — with integrally quan- 
tized magnetic flux ^ (n = 0,±1,±2 • ••; vorticity )0, 
hence the vortex could be anyon. The starting Ginzburg- 
Landau Lagrangian with the Chern-Simons term is an 
extended version of the Chern-Simons-Higgs model, so 
to say, the non-rclativistic Maxwell-Chern-Simons- vector 
Higgs model. We also discuss electromagnetic property 
of the vortex. In the region r < A (r is the distance from 
the vortex core, and A is London penetration depth), 
the vortex has an electric field from the Chern-Simons 
term and it would be expected that non-trivial electro- 
magnetic phenomena occur. Sr2RuOi is known as a 
type II superconductor J23 therefore it has a stable vortex. 
It is interesting if these exotic properties of the vortex 
are observed experimentally in the superconductivity of 
Sr2Ru04. 

The paper is organized in the following manner. In sec- 
tion II, we review the derivation of the Ginzburg-Landau 
Lagrangian with the Chern-Simons term in chiral p-wave 
superconductor. In section HI, we study charge and an- 
gular momentum of a vortex. In section IV, we discuss 
electromagnetic property of the vortex and find that an 
electric field is induced near the vortex core. Summary is 
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given in section V. We use the natural unit {Ti = c 
in the present paper. 



II. DERIVATION OF THE CHERN-SIMONS 
TERM AND THE GINZBURG-LANDAU 
LAGRANGIAN OF CHIRAL P-WAVE 
SUPERCONDUCTOR 

Following RefEi, we review how the Chern-Simons 
term is induced in the Ginzburg-Landau Lagrangian of 
chiral p-wave superconductor such as Sr2RuOi in this 
section. At first, we consider one superconducting layer, 
and we extend our discussion in multi-layer system. 

Fermionic Lagrangian of superconductor is written as 

1 



.^t I e(P + eA) A(p) 



At(p) -e(p-eA) 



(1) 



where 



is Fermion field in Nambu- 



Bogoliubov representationcJ with "isospin" a, and is 
the usual representation for Fermion field with spin index 
s- {T)ap is Pauli matrix with isospin indices, and e(p) 
is the kinetic energy of Fermion and is measured from 
Fermi energy ep. For simplicity-. we choose 1 band model 
with cylindrical Fermi surface,cl therefore, e(p) is writ- 

2 

ten as e(p) — ep with electron effective mass mg. 

A(p) ~ Vpair < V'V' > is the gap function (<•••> de- 
notes the expectation value of operator in the ground 
state and Vpair shows electron-electron interaction by 
which pairing state arises) and it is written as 

A(p) = iaa2 ■ d{p), (2) 

in the p-wave (spin-triplet) pairing statcEl (fss' is Pauli's 
spin matrix and d is so-called c?-vector. In the supercoa- 
ductivity of Sr2RuOi, proposed d- vector is written asEl 



d(p) 
?7(x) 



ZT] ■ p, 

fe(x),ryy(x)), 



(3) 



where 77 is a vector order parameter field and it has an 
expectation value written as 



< V >= hoKl,*), 



(4) 



which shows chiral p-wave pairing state. For a while, we 
neglect the spatial dependence of 77 and put 77 =< 77 >. 
Gauge transformation of ^' and Ap is written as 



A,. 



A,. 



(5) 



For convenience, we take doAo — c^diAi = gauge, where 
Cs is the velocity of U{1) Goldstone mode, and fix the re- 
dundant gauge degree of freedom which corresponds to 
the degree of freedom of U{1) Goldstone mode. 



It is known that the massive particle gives no higher or- 
der corrections to the Chern-Simops term, but the mass- 
less particle may give a correctionBa ICp find out a higher 
order correction, the authors of Rcf.E3 considered a cor- 
rection to the induced Chern-Simons term by U{1) Gold- 
stone mode. It was shown that U{1) Goldstone mode 
does not give a correction to the induced Chern-Simons 
term and its degree of freedom corresponds to the redun- 
dant gauge degree of freedom of the system. 

From Eq.(^, we can obtain Fermion propagator and 
calculate Fermion loop. Before doing that, to simplify 
notations, we introduce a vector ^(p) as 

/ ReA(p) \ / \r,o\p 
5(P) = -lmA(p) = -\vo\Pyia3a2 . (6) 

V 6(p) ; V e(p) J 

By using g{p), Eq. (|l|) is rewritten as 

C = *t[^5o-T•.g(p)]vE'-j^A^ 
jo(q) - e^*^3*p+q, (7) 

j(q) = j^^Hq)+j^'Hq), 

where j'^^'^ is the paramagnetic current and j is the 
diamagnetic current. They are written as 



j(P)(q) = eE*Ui53(p+f)}* 



p+qj 



j^'Hq) = ^E*^A(q)r3Vl/p+q. 



(8) 



We regard the interaction jf^A^^ as a perturbation and cal- 
culate correlation functions. Fermion propagator Gp{t — 
t') =< T^p{t)^l{t') > (T means the T product) ob- 
tained from Lagrangian Eq. (0) is written after the 
Fourier transformation with (2-t-l)-dimensional space- 
time momentum = (po, p) as 

G{p) = (po - T ■ gip))-' 

where iS stands for a infinitesimal imaginary constant 
and it fixes the boundary condition for the propagator. 

The Fourier transforms of density-density, density- 
current, and current-current correlation functions in the 
lowest order of the perturbative expansion are written as 

7roo(9) = -i{F.T.) < Tj^{x)jo{x') > 
d'^p 



Po 

pI- 



.Tr [T3G(p + g)T3G(p)] 



(27r)3 

7ro.(g) = -i{F.T.) < Tjo(a;)jf ^ (x') > 
d'^p 



(10) 



-le 



T,G{p + q){-^)G{p) 



dp^ 



(11) 



-t{F.T.)<Tj^P\x)&\x')> 



-le 



d^P gsjp) .93 (p) 
(27r)3 dp' dp> 



Tr [G{p + q)G{p)] 



= Oiq'), 



(12) 
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and a contribution from the diamagnetic current in the 
same order is 



me 



(13) 



where x denotes the coordinate in (2+l)-dimensional 
spacetime, {F.T.){- ■ ■) means that we make the Fourier 
transformation of following terms, Tr means the trace 
about isospin and real spin indices, and pe = ^"^'^ is 
Fermion number density. Since we are interested in low 
energy and long distance region, we neglect 0{q^) terms. 
The constant axy will be the coefficient of the Chern- 
Simons term, and it is written as 



'xy 



^2 



(Pp tr[g- {dg x dg) - gz{dg x dg)^ 



i2ny 



tr[y-g\-^ 



(14) 



here tr means trace about real spin indices and 9 = 
The first term in Eq.(14) is a topological im/aEiant. 
Volovik argued this type of topological invariant .llZlcj The 
second term is not a topological invariant. The second 
term appears bpsause U{1) gauge symmetry is sponta- 
neously broken .Ej 

By substituting Eq. (|^) into Eq. (|l^), we obtain the 
value 



e 

4^' 



(15) 



which coincides with the fine structure constant. In this 
case, the non-topological term in Eq. ( p^ vanishes ac- 
cidentally. In general, the non-topological term does not 
vanish. Actually, if we calculate (Jxy in the tight bind- 
ing scheme, the non-topological term has negligibly small 
value axidjixy is approximately equal to the fine structure 
constant .EEl 

From Eqs. (^, (0), (|l2|), and (|l|), the low energy 
effective Lagrangian in the case that the vector order pa- 
rameter field is constant (i.e., 77(x) =< rj >= |77o|(l,i)) 
is written as 



eff. 



-iw,. + i(^A^-m^A2) 



(16) 



The second term in the r.h.s. of Eq.([16D contributes 
Thomas-Fermi screening, and the third term contributes 
Meissner effect. The parameters m and Cg are determined 



/Pee^\l/2 



(^) 



1/2 



(22_)i/2 where 



Pe and v-p are mass, number density and Fermi veloc- 
ity of electrons, respectively. A = is the penetra- 
tion depth for the magnetic field in the usual case {i.e. 



0). The combination of the forth term and the 



fifth term in the r.h.s. of Eq.(16) is the Chern-Simons 



term in OqAq — c^diAi — gauge. 

Next we take into account spatial dependence of or- 
der parameter fields. Actual superconductor is spatially 
3-dimensional object which consists of layer of the super- 
conducting plane. Therefore, we extend above discussion 
to layered system. In low energy and long distance re- 
gion, it is sufficient to consider lower mass dimensional 
terms. Except for the Chern-Simons term, the quite gen- 
eral phenomenological Ginzburg-Landau free energy of 
the chiral p-wave superconductor with tetragonal symme- 
try which corresponds to the layered perovskite structure 
of the_crystal lattice has been proposed by Sigrist and 
Ueda.L3 We combine both of them. Lagrangian would 
have the following form ; 

+\DoV\'' - Ki{\Dxm\^ + \Dym\'') 

~K3{{DxmnDyV2\) + c.c.} 

-Ki {{Dxri2T{Dyjn\) + c.c.} - Xsl^zr/P 

-A\T^\^-PM''~P2W -^v? 

-Iiz\rii?\ri2? + C, (17) 

where = + 2ieA^, rip is a number of the super- 
conducting plane per unit length along the z axis (the c 
axis), and A, Ki, and are phenomenological constants 
with appropriate dimensions. 

To simplify the calculations, we use 



K^ = K2 = ci, = -Ki 
A = f + 6, /32 = -6, /33 = 0, 



A 



(18) 



, (0<a,6). 



These assignments of the piienomenological parameters 
in Eq. ( p^ ) seems artificial, however, charge and an- 
gular momentum of a vortex do not depend on these 
phenomenological parameters as will become clear in the 
next section. Eq. ( [T7| ) reduces its form with cylindrical 
symmetry to 

Cg-l = np^e'^^iA^d^A, + A^d^A^) ~ -^F^^^F^, 

-2i\v\'-\Vo\'r-b{\v\'-W^v\'), (19) 



where D 



// 



y). The Ginzburg-Landau La- 
(|19[) is an extended version of the 



grangian Eq 

Chern-Simons-Higgs model, so to say, a non-relativistic 
Maxwell-Chern-Simons-vector Higgs model. 
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III. VORTEX WITH FRACTIONAL CHARGE 
AND FRACTIONAL ANGULAR MOMENTUM 
IN CHIRAL P-WAVE SUPERCONDUCTOR 

In this section, we study charge and angular momen- 
tum of a vortex in chiral p-wave superconductor. As we 
mentioned before, charge and angular momentum take 
fractional values. From these facts the vortex in the 
present case could be anyon.Ej 

Suppose the magnetic field is directed to the z-axis. We 
consider static and cylinder-symmetric solution around 
the z axis and we neglect time and 2;-dependence. Eqs. 
of motion in V • A = gauge are written as 

V • E = -Se^lrypAo + Upa^yB, 
(V X B), = -Se^c'M^Ai + 2iecl{ri*d^rj - r^a.r;*) 

+2iec1-f€ij3dj{r]* x 77) + npa^yeijidjAo, 

(z = l,2) 



+ ^i\v\'-\rio\')Vk 



(20) 



where E* = and B = ~F^^. We see that there are 
two typical length scales in these equations. One of them 
is A = (Sc^c^It^oP)"^^^ for gauge fields and another is 
^ = Cs{y^\r]o\)^^ for T]. X coincides with so-called Lon- 
don penetration depth of the magnetic field in the case 
o'xy = 0. We call A also "penetration depth" in this 
paper. ^ is the coherent length. The Ginzburg-Landau 
parameter k is defined as k = A/^. It distinguishes type 
I superconductor (k < -^) and type II superconductor 

(-^ < k). Type II superconductor has a stable vortex 
solution but type I does not. The chiral p-wave supercon- 
ductor Sr2RuOi has values A ~ I8OOA and k = 1.2 for 
the magnetic field directed to the c axis OiLjthe crystal, 
which corresponds to z axis in this paper.t3 Therefore, 
Sr2RuOi is a type II superconductor and it has a stable 
vortex. 

We use cylindrical coordinate and consider an Ansatz 
written as, 



77 = \sA^p^^r),zp2{r))e^ 
^0 = 25Aao(0> 



(21) 



A. 



2er 



{n~ae{r)), 



which shows a situation that a vortex with vorticity n is 
located at the origin. For the total energy of the system 
to be finite, boundary condition at infinity is written as 



Pi{oo) = 1 
ao(oo) = 0, 
a0{oo) = 0. 



i = 1,2, 



(22) 




0, 
const. 

n. 



(23) 



are required to exclude singularities of the fields. 

First, we consider the asymptotic behavior of Ansatz 
Eq.(|l]). As we mentioned before, we see from Eq.(|20|) 
that a typical length scale of is ^. Hence we regard 
Pi ~ 1 when r is sufficientbi, larger than ^ and behavior 
of flo and ae are written aaij 

ao{r) ~ C+Ko( — ) + C-np(j^y\clKo{-), 

ae{r) ~ -j |c+c,npa,yAKi(-^) - C_Ki(^)| , (24) 

where K/ (r) is the ^-th order modified Bessel function and 
(7+ and C_ are numerical factors. The power behavior 
of the fields near the origin is written as 

ao ^ ao + !^4^)(r/A)2 + ©((r/A)^), 

ag^n + af\r/Xf+Oiir/\)^), 
P.'-pl"^(^/A)l"l+0((r/A)l«l+2), (25) 

where ao, a''g\ p^"'' are numerical constants. 

By using asymptotic behavior of the fields, we calculate 
the charge and the angular momentum of the vortex. In 
the rest of this section, we consider pure 2-dimensional 
space to avoid complexity in the arguments. It corre- 
sponds to putting rip = 1 in our calculation. 

The vortex charge is obtained by integrating charge 
density of the matter, such as 

d?x {V • E - a^yB} 



-2^a,,(rA,)|™ 
2n7r ne 



(26) 



Boundary conditions at the origin 



by using Eqs. (|2^), (|l|), (|2|) and (|5|). We can see 
the flux is attached to the charge, i.e. the vortex 
is an object such as flux-charge composite. Volovik has 
mentioned the vortex charge caused by the Chern-Simons 
term in superfluid '^He-A film. His result has half value 
comparing with our result Eq. ( p6| ) because the factor 
2 was missing when thcj-vprtex charge was derived from 
the Chern-Simons term.Ea 

It was argued by many authors that the change in the 
density of electrons in the vortex core due to the spatial 
dependence ot-the orderparameter also gives the charge 
of the vortexc3 and we call it "the regular charge" . Re- 
cently, Volovik pointed out that the occupation of the 
zero-energy bound states of electrons, which existin the 
vortex core of the chiral p-wave superconductor could 
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contributes the charge of the vortex.iil@ We call it "zero- 
mode charge" . It is important how to distinguish these 
two charges and the fractional charge Eq. ( p6| ) comes 
from the Chern-Simons term. It can be done, even in the 
case that the regular charge or the zero-mode charge are 
of order e, by comparing the vortex charges of vorticities 
n — 1 and n — —1 because the fractional charge Eq. ( |2^ ) 
changes its sign but other two charges do not depend on 
the sign of vorticity. 

The definition of the vortex angular momentum is 



J 



d^xr X v. 



(27) 



V is the momentum density, which is the generator of the 
translation. It is defined to be gauge invariant such as 



dC 



-D,ri + {h.c.) + 



[DoriY D^T^ + Dar]{D,ri) 



e^EjE. 



(28) 



By using Eqs.(pO|), ( p^ and (p5|), angular momen- 

tum Eq.(p7|) is calculated as follows; 

J = cfxrpg 



I zer J 



-2vr^ dr-dr{n^ae{r)Y 



n 

16' 



(29) 



We see the vortex in our system Eq. (|l9|) has fractional 
charge and fractional angular momeptum same as the 
vortex in Chern-Simons-Higgs model.E^I Since the charge 
and the angular momentum are proportional to a^y, we 
see that the Chern-Simons term creates these non-zero 
fractional values. These values are topological and de- 
pends only on asymptotic behavior of the gauge fields, 
although there is Maxwell term —jF^^F^^, and the Higgs 
field is vector-like. 

The results Eq. ( p6| ) and Eq. (|29| ) would be general 
and would not depend on the form of Ansatz and phc- 
nomenological parameters. Actually, in amconventional 
superconductors there is another AnsatzLd which is dif- 
ferent from Eq. (|2l|), and there are plausible phenomenct 
logical parameters for Sr2RuOi studied by AgterbergcZI 
and they are different from Eq. ([l^). The results Eq. 
( p6| ) and Eq. (^9|) would not be changed even if we use 
such a different Ansatz and phenomenological parame- 
ters. The reasons are written as follows; 

1. The definition of charge density and angular mo- 
mentum density do not depend on the phenomeno- 
logical parameters explicitly and always become to- 
tal derivative, since 



(a) charge density of matter j^^o-ttf^^ in Eq. 
(Eq) always satisfies the Maxwell equation 
^yatter ^ y • E - a^yB, and 

(b) the momentum density Vi defined by Eq. ( |2^ ) 
depends on the form of time derivative terms 
in Lagrangian Eq. (|9|) and these terms do not 
contain the phenomenological parameters. 

2. Asymptotic behavior of the fields is investigated 
without the Chern-Simons_term in Ref.Eil by us- 
ing Ansatz written in Ref.E3 and|-the phenomeno- 
logical parameters written in Ref.EZl. The essential 
property of the asymptotic behavior is the same 
with Eq. ( p^ ) and Eq. (p5[), and it would not be 
changed even if the Chern-Simons term is taken 
into account. 



Therefore, the results Eq. (|2^) and Eq. (|29|) would 
be valid for realistic systems such as superconducting 
Sr2RuOi. Evidence that the vortex we consider is a 
candidate of anyon can be seen in— the fact that vortex 
has fractional angular momentum£3 It is expected that 
transmutation of the statistics occurs. The further dis- 
cussion is needed to see the fractional statistics of these 
vortices experimentally. 



IV. ELECTROMAGNETIC PROPERTY OF 
VORTEX WITH FRACTIONAL QUANTUM 
NUMBERS 

In this section, we show numerical solution of Eq. ( |20| ) 
by using Ansatz Eq. (^ij) and discuss electromagnetic 
property of the vortex. We consider layered system again. 
In this case, Up ~ 10~^A~^. For simplicity, we assume 7 
in Eq. ([l^) is zero. In this case, there is an exchange sym- 
metry such as pi ^ p2 in the equation which is obtained 
by substituting Eq. ( pl| ) into Eq. (pO|). We introduce a 
real field 6 written as 



(r) = pi(r) = p2ir). 



(30) 



The magnetic field, the electric field and (f) with parame- 
ters K = 1.2 and Cs ~ (^)^/^ ^ lO^'^ in the case n = 1 
are shown in FIG. 1. Contrary to the usual vortex {i.e. 
the vortex in the case axy = 0), we found that an electric 
field in the radial direction is induced in the region ?- < A 
and its magnitude is about 1 Volt/meter. 

Let AE stands for energy difference between our vor- 
tex E(axy = 1^) and the usual one E{(Jxy = 0). There 
is a relation between them written as 



AE 



E{axy 



in I 



E{axy = 0) E{axy = 0) 



- 1 



(31) 



Therefore, energy increase caused by the existence of the 
Chern-Simons term is extremely small and the vortex 
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with fractional quantum numbers could exist in realistic 
systems. The magnitude of the vortex energy and also 
that of the electric field depend on the form of Ansatz 
and phenomenological parameters in Eq. (IT^), however, 
these calculation would have validity as an order estima- 
tion. It is interesting if some phenomena caused by the 
induced electric field are detected by experiments. 



V. CONCLUSION 

We have discussed a vortex in chiral p-wave super- 
conductor such as Sr2RuOi, where U{\), P- and T- 
symmetry are broken simultaneously. We have inves- 
tigated the vortex based on the Ginzburg-Landau La- 
grangian. The Ginzburg-Landau Lagrapgian of the sys- 
tem contains the Chern-Simons termjlj and it corre- 
sponds to an extended model of the Chern-Simons-Higgs 
model. We have found that the vortex in chiral p-wave 
superconductor has fractional charge — ^ and fractional 

angular momentum — -jg with integrally quantized mag- 
netic flux . These values are topologically stable and 
do not depend on the form of Ansatz and the phenomeno- 
logical parameters in the Ginzburg-Landau Lagrangian. 
Following the discussion in Ref.EJ, the vortex could obey 
the fractional statistics. We have also investigated the 
electromagnetic property of the vortex. We have found 
that the electric field is induced by the Chern-Simons 
term near the vortex core. Energy increase caused by 
the existence of the Chern-Simons term is so small that 
the vortex with fractional quantum numbers could ex- 
ist in realistic systems. It is interesting if these exotic 
feature is observed experimentally. 

Our discussion is valid for other chiral superconduc- 
tors, such as "anisotropic chiral p-wave" superconduc- 
tor in which the symmetry of the, wave function of 
Cooper pair is {smpx -l-«iSinpj,)-waveEJ and dx2-y2 +idxy- 
wave superconductorj23 since it has been shown that 
the Chern-Simons term is also induced in thej-Ginzburg- 
Landau Lagrangian of these superconductors. £3 
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FIGURE CAPTION 

FIG. 1. r dependence of (a) the orderpapameter ^, (b) the magnetic field directed to the ^-axis, and (c) the radial 
electric field. 
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